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Abstract 

We investigate the implication of unparticle exchange for the possible Sommerfeld enhancement 
in dark matter annihilation process. Assuming the unparticle exchange during WIMP collision, 
we solve the Schrodinger equation for the effective potential, and find that the Sommerfeld en- 
hancement factor is dictated by the scale dimension of unparticle as l/v 3 ~ 2du . Numerically the 
Sommerfeld enhancement could be O(10 — 10 3 ). 
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[Introduction] It has been known that our universe is made of not only the stuff of the 



r o, but also dark matter and dark energy 
11. Therefore, it must be one of the most 



standard model (SM) with the occupancy of mere 
with the abundance of 22% and 74%, respectively 
important issue to understand what the dark matter is and how to explore it by various 
observations. Hopefully, through high energy colliders such as the Large Hadron Collider 
(LHC) at CERN, we may directly observe dark matter soon. In the mean time we may also 
have the chance to probe dark matter indirectly by the study of the high energy cosmic-ray. 

Recently, the collaborations of PAMELA Q, ATIC Q, FERMI-LAT [l], HESS |s] and 
etc. have published quite astonished events in cosmic-ray measurements, in which PAMELA 
observes the excess in the positron flux ratio over 10 GeV till PAMELA'S observational limit 
of about 100 GeV, whereas others measure consistent anomalies in the electron+positron 
flux in the 300 — 1000 GeV range. Inspired by the new founds at the satellite, balloon and 
ground-based experiments, the excess could be readily ascribed to dark matter annihilation, 
even though there still are the possibilities of existing new young pulsars Although dark 
matter decays could be the origin of such anomaly, however, for making the lifetime as long 



as O(10 25 ) seconds, the extreme fine-tuning on the coupling of interaction 



ion [z, 8] 



cannot be 



avoided. For escaping the fine-tuning problem, hereafter, we will focus on the mechanism of 
dark matter annihilation only. In addition, the candidate of dark matter in our following 
analysis is regarded as weakly interacting massive particle (WIMP). 

Although the WIMP annihilation could be the source for the excess of cosmic-ray, how- 
ever, due to low reacting rate in the annihilating process, an enhanced boost factor of a few 
orders of magnitude, e.g. Sommerfeld enhancement [9 16] , has to appear during the anni- 
hilation of WIMP. Therefore, a new force carrier in the dark matter annihilation to dictate 
the enhancement is required. As an example, an interesting mechanism for the Sommerfeld 
enhancement is arisen from the light boson exchange between dark matter Kjl [ill , where 
the resulted interaction is Yukawa potential and the force carrier has the significant influ- 
ence in the range of Compton wavelength, denoted by aM x with a being the fine structure 
constant of the interaction and of order 10~ 2 . 

In this Letter, we study another kind of new force that may be alive in an invisible 
sector and dictated by scale invariant. As known that an exact scale invariant stuff cannot 
have a definite mass unless it is zero, therefore for distinguishing from the conventional 
particles, Georgi named the stuff as unparticle 18, 3]. Interestingly, it is found that the 
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unparticle with the scaling dimension du behaves like a non-integral number du of invisible 



particles 



Further implications of the unparticle to colliders and low energy physics 



could be referred to Refs. 



20 



22|. In order to concentrate on the Sommerfeld effect, here 



e, e.g. in our analysis, we 



23j. Although WIMP and 



we don't study the general effective interactions with unpartic 
have suppressed the interactions between unparticle and Higgs 
unparticle both weakly couple to the SM particles, however, there is no any reason to limit 
that the interactions between them should be weak. Therefore, when WIMPs collide each 
other with small speed, we speculate that the Sommerfeld enhancement could be arisen from 
the unparticle exchange during the collision, sketched in Fig. [TJ where the \ an d U denotes 
the WIMP and unparticle, respectively. And represents the (generic light) particle that 
might weakly decay to SM particles and the constraints on the couplings will be controlled 
by current observed fluxes of cosmic rays such as electrons, positrons, antiproton and etc. 
Here, its appearance is responsible for the possible connection between dark and visible 
sectors, but not for the Sommerfeld enhancement. Hence, we don't further discuss the 
detailed couplings to the SM stuff and the related issue for this decay. Our motivation is to 
understand that if there exists unparticle in the invisible sector, what are its unique character 
on the Sommerfeld factor and the differences from Coulomb and Yukawa interactions? 



FIG. 1: Dark matter annihilation with unparticle-mediated Sommerfeld enhancement. 



The paper is organized as follows: First, we derive the static unparticle potential, solve 
the associated radial Schrodinger equation with suitable boundary condition, and find the 
resultant formula for s-wave Sommerfeld enhancement. Then, we do the numerical analysis 
on the resulted Sommerfeld factor and present it as a function of involved parameters in 
two-dimensional contour plots. Finally we give the conclusion. 



[ Unparticle potential and its Sommerfeld factor] Although the nonperturbative Sommerfeld 
effect can be calculated by the combined contributions of a set of ladder diagrams shown 
as in Fig. [H however in the non-relativistic limit, the effect could be equivalent to solving 
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the Schrodinger equation with an effective potential, which is arisen from the single particle 
exchange. In the considered mechanism, here the exchanged particle is the unparticle. 



Following the scheme proposed in Ref . 18j , the interaction of WIMP to unparticle can be 
written as 



T X [c) xOu (1) 



where we have assumed that the WIMP is a Dirac (or Majorana) fermion, A is dimensionless 
parameter and Ay denotes the living scale of unparticle. For displaying the character of 
scale invariant stuff, we concentrate only on the scalar unparticle. Please note that when 
unparticle is realized in the framework of conformal field theories, the propagators for vector 
and tensor unparticles should be modified appropriately to preserve the conformal symmetry 



241. To obtain the unparticle potential in non-relativistic limit, we use the propagator of 

nn 



the scalar unparticle operator given by 

mm 

/e^<0|™ u (,) O „(0)|0> ^ i5 -^__L_, (2) 
where A du 



'2sin4/7r (-g2 
16tt 5 / 2 Y(d u + 1/2) 



(27r) 2 ^ T(d u - l)r(2d M ) ■ 

Combining Eqs. (CQ) and (jSJ), the four- fermion effective interacting term in momentum space 
could be expressed by 



x (c) x 



A \ 2 A 



Afj 1 - 1 J 2 sin d u ir {-q^f-du 



X {C) X- (3) 



HA 

By Fourier transformation and with q° = 0, the static unparticle potential in WIMP inter- 
action resulted by Eq. (j3J) is found by 

V(r) = -5 (4) 



with t = 2d 



u 



a 



( A 



2lT 2d U V \ d U~ l 



T(d u + l/2)T(d u - 1/2) ... 
fr = W • (5) 

Interestingly, we see that the power of unparticle potential is associated with the scaling 
dimension du and not an integer. 
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Since the unparticle potential is independent of the polar and azimuth angles in 
spherical coordinate system, the relevant piece for the Sommerfeld factor is the radial 
Schrodinger equation, read by 



d 2 ( £{£+!) 2/ia 2 



uu{r) = , (6) 



where we have used the nature unit with h = c = 1, E = k 2 /2fi with /x being the reduced 
mass of system and Rki{r) = u^(r)/r. Once we solve the differential equation with the 
proper boundary condition, the Sommerfeld effect associated with each angular momentum 
£ is obtained by 



S, 



(2£+l)\\d i R ke {r) 



r=0 



(7) 



\k\*£\ dr e 

Due to the power t of unparticle potential in r being not an integer, there is no hope to find 
a general close form for the solution. Nevertheless, since the required Sommerfeld factor is 
estimated at r = 0, our strategy for finding the solution is to look for a good approximation 
to extract the behavior of radial wave function at r ~ 0. 

Before discussing the solution to the differential equation, first we analyze the limit on 
t = 2du — 1. To avoid the crossing of the branching cut at 1/ sino?^7r in Eq. (J2D, we require 
1 < du < 2 (1 < t < 3). In order to further understand whether the upper limit of t can 
be bounded by the boundary condition when solving differential equation, we examine the 
case with t = 2, in which the potential could be exactly solved at r — > 0. Hence, by taking 
t — 2 and Uke ~ r a and considering r — > 0, from Eq. © we get 



a= l - + ^{£+ l -f-2 m . (8) 

Thus, the corresponding solution for Uke{r) can be expressed by Uke(r) = r a e lhr f (r) , where 
the function of f(r) is controlled by 

d 2 f df 
r—^- + {2a + 2ikr) -f - + 2ikaf{r) = . (9) 
dr z dr 

Compare to the confluent hypergeometric differential equation, xy" + (c — x)y' — ay = 0, in 
which the solution is confluent hypergeometric function of the first kind y(x) = \Fi(a, c; x) , 
we see that the solution to Eq. fl9]) can be found by taking a = a, c = 2a and x = —2ikr. 
For £ = 0, from Eq. ([8]) it is easy to find a < 1. Since 2a, —2ikr) — > 1 when r — > 0, 

as a result, the wave function Rke(r) is singular at origin. In other words, t = 2 leads to an 



ill-defined solution and should be taken as the new upper limit. Therefore, we find the new 
bound on the scale dimension of unparticle as 1 < du < 3/2 (1 < t < 2). 

After obtaining the allowed power of unparticle potential, next we are looking for a good 
approximation to get the proper solution for near origin. As used before, if we regard 
Ukt( r ) = r e+1 e lkr f(r) as a general solution to the radial Schrodinger equation, according to 
Eq. ([6]), the function /(r) should satisfy the differential equation 

r% + W + 1) + 2i*r)! + (2ik(£ + 1) + ^ f = . (10) 

Although this equation cannot be solved for arbitrary value of t, we find that if we concen- 
trate on small r, the terms associated with 2ikr and 2ik(£ + 1) can be safely neglected and 
the solution to the simplified differential equation can be found by 

f t {r) = d -(^ + D/ 2 (2-.) (2 _ i)^^ 1 )/^^^^^ 
,'2£-t + 3\ T [2i+1 2yfdr^l 2 

x r — - — ■ — j 



t J I 2-t ' 2-t 

where J u (x) = J(v,x) is the Bessel function and the asymptotic form for < x \Jv + 1 
is J u — > (x/2) v /T(v + 1). Hence, the solution to radial Schrodinger equation in the region 
of small r can be expressed by 

Ru{r) = ^ = rV fe 7t(r) , (11) 

where a suitable normalization for the wave function has to be further dealt with. 

Unlike conventional approach that the normalization of wave function for scattering pro- 
cess is chosen at r — > oo, in order to get a suitable normalization for the wave function at 
near origin, we have to consider an alternative method which could provide correct value 
of wave function at r = 0. We find that the purpose can be achieved by using integral 
form instead of directly solving the Eq. (Q. In other words, the solution for Uki{r) could be 
expressed by 



Uki(r) oc j t (kr)+ / dr'G{r,r')U(r')u M (r') (12) 
J o 

with U(r) = 2fiV(r) and 

, ; fo rr </, 

G{r,r) = { (13) 
n i{kr)ji{kr') jk for r > r' , 
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where G(r, r') is the Green's function of the differential equation. After compared to the 
partial wave expansion given by exp(ikz) = X/£o* + ^)je(kr)Pe(cos9)/kr, we find that 
the solution for r — > is 

uu = 3J ^P- + 3j ^r- J" dr'n e (kr')Mkr')U(r') . (14) 

Now we match both solutions at r pa by solving differential equation and the integration 
with Green function. It has been known that the Sommerfeld factor of Coulomb potential 
(i.e. t = 1) for s-wave can be obtained exactly as 

So = 2 T /V • (15) 

1 — exp(— 2aix /v) 

For the case of v <C 1, we have Sq ~ 2ira/v. Based on the known result, we find that in 
order to get correct approximation for Coulomb potential, matching condition for Eqs. ( ITT]) 
and (IT4"1) has to be adopted as 



Nir 



2^,^+1 



4 r°° 
= ]pM kr ) J dr'n e (kr')ji(kr)U(r') . 



Using the asymptotic results f(r — > 0) — > 1 and je(kr — > 0) — > ye+iy. (kr) e+1 , the normal- 
ization constant is found by 

= ^j0Yy kl+t ~ 2x " (16) 

with Xi = J °° dzji(z)n£(z) / z l . According to Eq. (J7J), the Sommerfeld factor for s-wave is 
obtained by Sq = 8a\Xo\fi t ~ 1 /v 2 ~ t , where we have used k = \iv. One can see that for t — 1, 
X = — 7r/4 and So = 2ira/v. Hence, the s-wave Sommerfeld factor induced by unparticle 
exchange is given by 

2A\ 2 , (fiY^-V \X 



s °-^l MvJ (17) 



[Numerical analysis and discussions] We now analyze the formula for Sq numerically. 
From Eq. (TTTj) . we see that the involved free parameters are A, dy, Ay and v. Unlike the 
case of Coulomb or Yukawa potential, the velocity-dependent factor appears by \/v z ~ 2du 
with 1 < du < 3/2. The Sommerfeld factor is increasing when du approaches to unity. 
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FIG. 2: Contour for Sommerfeld factor induced by unparticle-mediated as a function of A and du 
with v = 10~ 3 . The numbers in the plot denote the values of Sq. 

Consequently, if we focus on the maximum of So, then we find that So is insensitive to the 
values of /i and Ay, where they show up by (/^/A^) 2 ^ w_1 \ Note that to explain the excess 
of cosmic rays from both PAMELA and ATIC/Fermi-LAT, it has been known that the mass 
of dark matter is of order of TeV. On the other hand, in order to produce the unusual stuff 
at LHC, the interesting scale to form unparticle should be also at TeV scale. Therefore, 
without lose of generality we set \x = Ay in our numerical estimates. 

In order to see the influence of remaining parameters, we will fix one parameter in turn 
when we make two-dimensional contours as a function of remaining two parameters. First, 
we analyze the contour as a function of A and du when the speed of WIMP is fixed. Although 
dark stuff weakly couples to visible particle, however, the interaction between invisible par- 
ticles may not be small. Accordingly, we set the allowed range for A be within one order of 
magnitude, i.e. 1 < A < 10. As a result, the contour for Sommerfeld factor So induced by 
unparticle contributions as a function of coupling A and scale dimension du with v = 10~ 3 is 
shown in Fig. [2j where the numbers in the plot stand for the values of Sq. By the figure, we 
see that So is increasing while du is decreasing. Furthermore, the value of du for S ~ (9(10 3 ) 
can be somewhat larger, when A is away from unity. It is clear that Sommerfeld factor could 
be as large as (9(10 3 ) by the new force mediated by unparticle. 

Secondly, we also study the behavior of So in A and v (in units of 10~ 3 ) by fixing the value 
of du- We present the resultant contour with du — 1.1 in Fig. [3] Clearly, with the value of 
du that is close to unity, Sq can easily reach O(10 3 ) even with A ~ 1 and v ~ 10~ 3 . On the 
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other hand, Sq of O(10 3 ) can be preserved when v and A both are increasing simultaneously. 

Finally, for understanding the dependence of speed of WIMP and scale dimension of 
unparticle on Sommerfeld factor, we checked the contour for So as a function of v (in units 
of 10 -3 ) and du with A = 1. And we found that Sommerfeld enhancement of (9(10 2 ) 
for explaining the excess of observed cosmic-ray through WIMP annihilation can still be 
achieved with v ~ C(10~ 3 ) and A = 1 while 1 < du < 1.1. 

As summary, inspired by the excess of electrons and/or positrons observed at PAMELA, 
ATIC, Fermi-LAT and etc., the issue of dark matter annihilation for the solution is revived 
and studied broadly. Besides a new mechanism is needed for the production of excessive 
cosmic-ray, usually we also need a new force that interacts between dark matter for over- 
coming the low cross section when dark matter collides. For studying the possible new force, 
we have investigated the impact of unparticle which is ruled by scale invariant. 

We find that the power of the static unparticle potential in r, resulted from the exchange 
of scalar unparticle, is non-integral number and it depends on the scale dimension of 
unparticle, expressed by \/r 2du ~ l . By the boundary condition for wave function at r = 0, 
the upper bound of scale dimension is found by c?^| max < 3/2. By looking for the suitable 
boundary condition and matching condition for the radial wave function at r — > 0, we 
are led to the Sommerfeld factor for s-wave collision, in which the result could return 
to that of Coulomb potential. However, unlike Coulomb or Yukawa potential, the speed 
dependence is related to scale dimension of unparticle and dictated by 1/ v 3 ~ 2du . Although 



FIG. 3: Legend is similar to Fig. [21 where du = 1.1 and the variables for the contour are A and v 
(in units of 10 -3 ). 
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the Sommerfeld factor is associated several parameters, however, our results are only 
sensitive to the parameters A, du and v. According to our numerical calculations, we 
conclude that with the allowed range of free parameters, the Sommerfeld enhancement 
induced by unparticle-mediated effects could be (9(10 — 10 3 ) and the factor could provide 
the necessity for enhancing the cross section of dark matter annihilation. 
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